The Effect of the Tilted Field in the Fractional Quantum Hall Systems: 
Numerical Studies for the Solid-liquid Transition 
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We construct a generalized Laughlin-liquid wave function and a variational electron solid wave 
function when the magnetic field is tilted. The energy of the liquid state is evaluated by Monte 
Carlo methods while the energy of the solid state is calculated by the optimization. Comparing 
these two energies for a given tilted angle 8, it is seen that the critical filling factor u c of 6 of the 
solid-liquid transition increases as the tilted angle. The implication to the experiment is that: i) the 
insulating phase may harder be melt as v — * 1/5 such that the width of the valley of the longitudinal 
resistance may become narrow as the filed is tilted; ii) it is expected that even in the vicinity of 
v — 1/3 for the electron system in the presence of the tilted field, the insulating phase may be 
observed. 



PACS numbers: 73. 20.Dx,73.40.Kp, 71.45. Nt, 

After the fractional quantum Hall effects (FQHE) 
was revealed, the measurement in a tilted field has be- 
come an established technique. Mostly, the tiled field is 
used to explore the spin related effects || . It also affected 
the transport of the bilayer system due to the in-plane 
field couples to the psuedospin degrees of freedom ||. 
The high Landau level stripe phase is strongly influenced 
by the tilted field [Q . For the spin polarized system, the 
tilted field effects are produced by either the finite thick- 
ness of the layer []5| or the subband Laudau level coupling 
||. For the later case, only a numerical diagonalization 
for a small system was carried out. One of focuses in 
the present work is trying to construct a Laughlin-type 
many-particle wave function when the magnetic field is 
tilted. 

In the early days that Laughlin's wave function was 
established, a set of the numerical works to check Laugh- 
lin's original wave function was presented 10-^]. Com- 
pared to the energy of the electron solid state, which is 
estimated by a variational method Ilfl] , it was found that 
the critical filling factor v c is around^r . Considering the 
Landau level mixing, it was said that the electron solid 
state may appear in a much larger filling factor pd| , ^2[ . 
The finite layer thickness can soften the short range be- 
havior of the Coulomb interaction and thus leads to favor 
the competing electron solid jl^,[Ll],[L4j . However, at the 
exact v = 1/m, the Laughlin-Jastrow correlation for the 
electron solid may leads to the melt of the electron solid 
|^5| while the energy of the liquid state at an arbitrary 
filling factor is higher than that from the interpolating 
between the odd-denominator filling factor |ll| , which 
may explain an reentrant insulating phase near v = 1/5 
for the electron jQ| and v = 1/3 for the hole [[ll| . 

In this work, we plan to study the competition be- 
tween the solid and liquid states as vary of the tilted field 
through the subband Landau level coupling. We perform 
Monte Carlo simulations for a generalized Laughlin state 
in the presence of the tilted field. As expected, the liq- 
uid state energy is lowered as increasing the tilted an- 



gle. Meanwhile, we use a variational wave function for 
the solid state in the tilted field and optimize the corre- 
sponding energy. The solid state energy is also lowered 
as increasing the tilted angle, whose magnitude lowered 
is always larger than that of liquid state energy. We 
can see that the critical filling factor for the solid-liquid 
transition may be considerably lifted. In this work, we 
do not consider the Laughlin-Jastrow correlation to the 
electron solid yet, which leads to an instability of the elec- 
tron solid exactly at v = 1/m for to an odd integer |l5| . 
Thus, although the tilted field effect may not destroy the 
liquid state at the exact filling factor 1/5, the insulator- 
quantum Hall liquid transition may become harder to be 
melt as v — > 1/5 such that the width of the valley of the 
longitudinal resistance may be narrowed as the tilted an- 
gle increases. Taking the energy lowered by the Laudau 
level mixing in the zero tilted field [O] as an estimation, 
one may even expect that the insulating phase will be 
observed near v — 1/3 for the electron systems in the 
presence of the tilted field. 

To consider the subband effect, we assume a quasi- two- 
dimensional electron gas confined by a harmonic poten- 
tial with V(z) — ^m e fl 2 z 2 , where m e is the electron band 
mass and f2 is the frequency of the harmonic potential in 
the perpendicular direction z. The applied magnetic field 
is tilted by an angle 9 to the two-dimensional x — y plane. 
The single-particle problem has been solved in literature 
Jl8| , [l9t . To construct a many-particle Laughlin-type state 
for the convenience of the numerical calculation, we take 
an alternative way. Working in the symmetric gauge, 
although the rotational symmetry of two dimensions is 
broken by the in-plane field, there is still a conserved 
quantity, = b^b^ no matter whether the magnetic field 



is tilted, where £ 



iy and — -1=(— <9? + |£) is the 
generator of the magnetic translation. The eigen state of 
thus is the same as that in B\\ = 0. The eigen state 
of the single-particle Hamiltonian for Bu ^ 0, then, can 
be the linear combination of the degenerate states with 
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a given eigen value mi, of L^, 
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where the unit of the magnetic length Ib ~ 



has been used and = 775 (— d$ + g£) and = 
^|(— d z i + z') are the harmonic operators with the fre- 
quencies lj c — and Q, respectively; z' — SV/ 2 z with 
Q = Q/u} c , Vl 2 = fl 2 + u> 2 and uj\\ == The combi- 

nation constants C™*™| can be determined by a standard 
method to find the common eigen states of two com- 
mutated operators. For the ground state, one denotes 
C°?n 2 = Cnm 2 , for example, C w = C i = and 



w c + £1 - w + - uj- 
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(2) 



for Coo is determined by the normalized condition. Other 
C nin2 can be deduced by recursive relations. Here 



c4 = ±(ft 2 +w 2 )± iyX^ 2 -u 2 ) 2 + AQ 2 fluj c are the cor- 
responding frequencies to the diagonalizcd harmonic op- 
erators. 

For the liquid state, following Laughlin's track in writ- 
ing down his wave function, the required incompress- 
ible wave function for the odd denominator filling factor 
v = 1/m has the form 

* ro (fi, ...,f 2 ) = IJ(E cw^r^D 

xIIte-y m e-^- (ly! ^ !) 



P f (ri, ...,fjv)|m), 



(3) 



where |m) is the common Laughlin state. To see the 
incompressibilty, one defines the electron density for this 
state 



p e (f) = (m\P^2s 3 (f - n)P^\m) I (m|PPt|? 



(4) 



Since H^E^f - n)\\m) = D r (rn\m) = (see §)) 
, it is easy to show p e (fO = (fn\ J^i S 2 (r — fi)\m) / (m\m) , 
which exactly equals to the density of the Laughlin state 
for B\\ = 0. Thus, the state is incompressible. 

In the same spirit, for the solid state, the variational 
wave function reads 

^ c {n) = P^ c {ii)e-^A\ (5) 
where ip c is the variational wave function fllCfl for Bn = 
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x exp 
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where Ui is the site fluctuations of the electrons away 
from the equilibrium sites Ri and Bij are the variational 
parameters. Similar to Lam and Girvin's, the optimal 
parameters, in the momentum space, are given by 



Bq = e 



W„ UL(q)-UJ T (q) 



where q is a phase factor and u>L{q) and u!r(q) are the 
effective longitudinal and transverse phone frequencies, 
which are determined by the dynamic matrix T> p0| , ^2| . 
Since the field is tilted, T> in the present case is related 
to an effective potential 
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n,2 — 0;ni ,1 — 1 



where M = E ni nA 



The energies we shall compare for the solid and liquid 

states are 
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l*c,m)- 



(9) 



Since the analytic single-particle wave function in the 
symmetric gauge is not ready yet, C nin2 has to be trun- 
cated [^0|. Fortunately, we find that for the practical 
value of Q/uJc and a modest tilted angle, the first order 
approximation, i.e., considering Cn-contribution only, 
may be a very good one. The expect energy of the ap- 
proximation single particle state compares to the exact 
ground state energy only a few percentage higher 
Thus, we consider the subband and the tilted field con- 
tribution from C'n only. 

For the liquid state, we evaluate the energy in (Q) by 
Monte Carlo methods. The detailed simulation technique 
has been described by Morf and Halperin || for the disk 
geometer. We consider a finite system. The positive 
charge neutralizing the electron is uniformly distributed 
on a disk of radius R m = \2mNlB- The energy calcu- 
lated, at Monte Carlo step s, is 
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where rf^ is the position of the z-th electron in the s-th 
Monte Carlo step; U(fi) is the energy of the i-th electron 
interacting to the neutralizing background: 
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where p e = -^W 1S ^e avera S e electron density. In the 
Monte Carlo simulation, we have used the unit of the ion- 
disk radius R = {-Kp^ 1 ! 2 . The probability distribution 
we used in the Monte Carlo simulation has included the 
tilted field effect, i.e., 



\* m \* = 



n(E^)4r4rn&-&) r 



x e 
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(12) 



In the real calculation, the perpendicular component z' 
can always integrated out by hands. We identify this 
energy of the s-th Monte Carlo step as the liquid state 
energy per electron. The simulations were done for the 
finite systems TV = 20,30,42 and 71. For N = 20 and 
30, this energy was computed by 50 independent Monte 
Carlo simulations, each of which consists of 170, 000 
Monte Carlo steps. For N — 42 and 71, the correspond- 
ing numbers are 30 and 110, 000. The results for E m are 
fitted by the polynomials 



E% = -0.782133i/ 1/2 + 0.148i/ 1241 - 0.009!/ 218 , (14) 
Ex. = -0.782133^ 1/2 + 0.144;/ 1 - 241 - 0.017i/ 218 . 

4 

For n/uj c = 0.4, the fitting results are 

Ex. = -0.782133^ 1/2 + 0.146i/ 1241 - O.OIO^ 218 , 

Ejl = -0.782133i/ 1/2 + 0.138jA 241 - 0.015zy 218 , (15) 

E« = -0.782133i/ 1/2 + O.llOzA 241 - 0.026z/ 218 . 

4 

The first one in ( |l4"| ) is corresponding to 9 — 0, which 
is consistent with the result of Levesque et al ||. One 
observes that the exponents in the fit formulas are not 
dependent on the tilted field; Neither are almost the co- 
efficients of v 2AS except in the last one of (15) which may 



be because the approximation single-particle wave func- 
tion is a little bit poor in this parameter (Ae/e ~ 13%). 

Now, we turn to the electron solid. The strategy we 
take is the same as Lam and Girvin's M] : Expanding the 
Coulomb interaction to the second order of the fluctua- 
tion Ui and determining the variational parameter B^. 
Then, calculate the energy (||) exactly. To the harmonic 
approximation, the harmonic term E c j t of the energy is 
given by 



E r 



E n 



A + Bx/VN + B 2 /N, 



(13) 



N ^ 

q 



(u) T {q) +uj L (q)f 



where and in the following, the energy unit is e 2 /els- 
We identify the zero-th order coefficient A as the result 
of thermodynamic limit. The values of A(v) of 9 for 
Q/lj c = 0.8 and 0.4 are listed in Table I. 



= ^E( Tr ^) + ^ 4Det2? (9) 



(16) 
(17) 



where T> is the effective dynamic matrix and q is the wave 
vector in the Brillouin zone. Comparing to the uncorre- 



Table 1(a) The thermodynamic limit of E m for Q/uj c = 0.8. lated energy E UCyh = ^ X^( w t(<?) + ^1 (<?)), the fluctu- 
ation is lowered the energy because the transverse modes 
are favored while the longitudinal modes are disfavored 

s 

The numerical results of Eg(v) are as follows. For 
Q,/lo c = 0.8, one has 
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A(l/3) 


A(l/5) 


A(l/7) 


A(l/9) 





-0.4099 


-0.3275 


-0.2809 


-0.2499 


tt/7 


-0.4131 


-0.3292 


-0.2820 


-0.2509 


tt/5 


-0.4146 


-0.3300 


-0.2825 


-0.2511 


tt/4 


-0.4163 


-0.3307 


-0.2830 


-0.2514 


Table 1(b) The thermodynamic limit of E m for Q./uo c 
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A(l/3) 


A(l/5) 


A(l/7) 


A(l/9) 


tt/8 


-0.4151 


-0.3303 


-0.2827 


-0.2512 


tt/5 


-0.4176 


-0.3315 


-0.2835 


-0.2518 


tt/4 


-0.4258 


-0.3356 


-0.2862 


-0.2537 



For 9 = 0, A is independent of Vl/uj c , e.g., for v = 1/3, 
we have A = —0.4099, which is consistent with results 
by Levesque et. al || and Morf and Halperin 

According to the results of the thermodynamic limit, 
one can figure out the energies dependent on the filling 
factor v and the tilted angle 9. For Q,/u) c = 0.8, one 
can fit the energy Eq{v) by the following interpolation 
formulas 

E Q = -0.782133^ 1/2 + 0.166iA 241 - 0.009i/ 218 , 
Eil = -0.782133z/ 1/2 + 0.154zA 241 - O.OlOi/ 218 , 



E Q = -0.782133^ 1/2 + 0.2430z/ 3/2 + 0.1610^ 5/2 , 
E f = -0.782133^ 1/2 + 0.2200i/ 3/2 + 0.1598^ 5/2 , 

E f = -0.782133^ 1/2 + 0.1971J/ 3 / 2 + 0.1596^ 5/2 , 

E* = -0.782133^ 1/2 + 0.1728i/ 3/2 + 0.1594i/ 5/2 . 

4 

For n/cj c = 0.4, it is 

Ez. = -0.782133^ 1/2 + 0.1926i/ 3/2 + 0.1595^ 5/2 , 
E% = -0.782133^ 1/2 + 0.1619i/ 3/2 + 0.1593^ 5/2 , 
Ejl = -0.782133^ 1/2 + 0.0992^ 3 / 2 + 0.1592;/ 5 / 2 , 



(18) 



(19) 



The first one in ( |l8|) corresponds to the energy of 9 — 0, 
which is consistent with Lam and Halperin's result Jl0| . 
We see that again the coefficients of v^l 2 are almost not 
affected by the tilted field. 

Using the fitting results flT|) and (||), (|l|) and |l|), 
one can obtain the estimated critical filling factors of 9 
(Table II). 
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Table 11(a) The critical filling factors for Q/uj c = 0.> 



9 


AE e (l/3) 


v c 





-0.015 


O.uo 


tt/7 


-0.013 


f! 4^ 


tt/5 


-0.010 


5.58 


tt/4 


-0.007 


4.54 


Table 11(b) The critical filling factors 


for Q./uj c = 0.4. 


9 


AE e {l/3) 


v c 


tt/8 


-0.009 


5.46 


tt/6 


-0.007 


4.45 


tt/4 


-0.002 


3.46 



In Table II, AEe(l/3) is the energy difference between 
the liquid and solid at u = 1/3, whose magnitude reduces 
considerably as the tilted angle increases from to 7r/4. 
In this work, we did not check the Laughlin-Jastrow cor- 
relation to the electron solid. It has been pointed out 
that the solid may be unstable at the exact v = 1/5 due 
to this correlation [jl5). Thus, our variational calculation 
may not be enough to completely understand the physics 
at the exact odd denominator filling factor 1/m. Due to 
the instability, the Hall plateau may not disappear. How- 
ever, our calculation indicates that the insulating phase 
may harder be melt as v — > 1/5 such that the width of the 
valley of the longitudinal resistance may become narrow 
as the tilted angle increases. 

On the other hand, although our result did not show 
the critical filling factor is shifted to near 1/3, the magni- 
tude of AEg(l/3) reduces remarkably as the field is tilted. 
According to Ref . |l2j , after accounting the Laudau level 
mixing, the energies of the solid and liquid states for 
9 = are lowered about 0.8% and 0.1% respectively at 
v ~ 1/3 and r s = ^ h '^ ° 2 ~ 2. The lowering of the en- 
ergy of the solid is larger than that of the liquid about 

2 

0.002^—, which can nearly overcome the energy differ- 
ence AEg(l/3) for large angles 9 in Table II. If one as- 
sumes the energy lowering due to the Laudau level mixing 
is independent of the tilted field and considers the energy 
of the liquid state for v being nearly but not exactly 1/ 3 
is in fact higher than our interpolating estimation p6[ , 
it is possible that at a large tilted angle, the energy of 
the solid becomes lower than that of liquid for v ~ 1/3. 
Thus it may be anticipated that in the vicinity of v = 1/3 
for the electron system in the presence of the tilted field, 
the insulating phase may be observed for a larger tilted 
angle. To verify this expected phenomenon, both further 
experimental and numerical works are required. 
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